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The present article represents part of the PhD. dissertation by C. Josser and§. We discuss 
the nucleation of quantized vortices in the nonlinear Schrodinger equation (NLS) for a flow 
around a disk in two spatial dimensions. It appears that the vortices are nucleated when 
the flow becomes locally (at the edge of the disk) supersonic. A detailed study of the phase 
equation for the complex field ip gives an Euler-Tricomi type equation for the stationary 
solutions below threshold. This equation is closely related to the one known in shock wave 
dynamics for gas. Then using solvability condition, we extract a time-dependent scenario 
for the evolution of the amplitude of the solution, which we, finally, relate to a known 
family solution of NLS which gives rise to a vortex nucleation. We also give a first order 
correction at the Landau velocity of nucleation, taking into account the geometry of the 
flow. 
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The dynamics of superfluids at zero temperature can be modelled by the Gross-Pitaevskii 
equation 0. This is a partial differential equation for a complex valued scalar field tp(x.,t) 
with dependence on the position in space, x and on time t. At zero temperature, there is 
no formal damping term in the equation; it is reversible in time (after complex conjugation) 
and has even a Lagrangian structure. Although many of properties are well-known and have 
been studied quite extensively, we shall first review some of them. This Gross-Pitaevskii 
equation can be seen as a nonlinear Schrodinger equation, and so shares many properties 
of the linear Schrodinger equation. With periodic boundary conditions, it has an uniform 
ground state minimizing the energy, the corresponding solution ip(x, t) depends periodically 
on time through a simple exponential factor. As the Gross-Pitaevskii equation is Galilean 
invariant, it is easy to construct a solution representing an uniform flow by boosting the 
rest state to a specified speed. With the same model, it is also possible to look at more 
complicated situations, as flows around obstacles. Two of the authors @JJ] have studied over 
the last years a 2D flow around a circular disk (among others), that is the solution of the 
Gross-Pitaevskii equation with a given uniform and constant flow speed and mass density at 



conclusions and further miscellaneous 
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I. INTRODUCTION 
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infinity, with a boundary condition on the surface of the disk. One striking observation has 
been that beyond a certain critical speed, the flow around the disk becomes time dependent, 
because vortices are emitted from the disk surface as shown on figure ([[]). 

In Ref. it was shown that the release of vortices from the boundary of the disk is 
a consequence of a transition from a locally subsonic to supersonic flow. In ordinary fluid 
mechanics, this would lead to the formation of a shock wave inside the so-called supersonic 
bubble. Nothing resembling this is possible in the Gross-Pitaevskii equation, because of the 
lack of built-in irreversibility, something that is necessary to balance nonlinearities inside 
the shock wave. As argued in Ref. ||, in the present model, the formation of shock waves 
is replaced -albeit in a rather loose sense- by the nucleation of vortices with a quantized 
circulation. In fact, these vortices are topological defects solutions of the Gross-Pitaevskii 
equation. In two space dimensions they are points where the complex field ip(~x,t) vanishes, 
making a ±2tt phase jump when one turns around the defect. These vortices are emitted 
when the local velocity becomes larger than the critical velocity making the flow locally 
supersonicQ 

The present work is devoted to study this transition. This is a rather intricate question, 
as anyone who looked at the paragraph §118, §119 & §120 of Landau and Lifschitz book 
on fluid mechanics || knows well: one has to solve locally the Euler-Tricomi equation for 
the velocity potential, that can be done in terms of the Airy function. In the present case 
there is a further complication with respect to Landau's analysis, which comes from the 
boundary conditions imposed on the disk, see §126 of Ref. 0. Let us outline the principles 
of our analysis; we are dealing with a problem, with two small parameters: the ratio of 
the microscopic intrinsic length to the disk radius, and the relative shift of the velocity 

1 The present problem has been studied in the one dimension case by Hakim [||] who has shown 
that the release of vortices is replaced by a periodic nucleation of one dimensional solitons. 
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near the pole of the disk to the critical value, a parameter called e.0 Far from the disk 
boundary, in terms of the microscopic distance, the usual hydrodynamic assumption holds: 
the quantum pressure term can be neglected, and one obtains an unique, nonlinear equation 
for the velocity potential^. This equation can be solved by expansion, by assuming that at 
the dominant order the surface of the disk is flat. In this limit, the hydrodynamic boundary 
conditions are satisfied with an uniform tangent velocity field. Corrections have to be added 
to this velocity field in order to take into account the curvature of the disk boundary. The 
first correction is trivial, and only the second one is crucial for the analysis. 

The second order corrections (for the velocity potential) is the solution of the Euler- 
Tricomi equation, with specific boundary condition. It happens that the small parameter 
e enters into this Euler-Tricomi equation plus the boundary conditions problem in such a 
way that it can be factored out by rescalings. Moreover, for e negative, that is for velocities 
slightly less than the critical speed, the velocity potential at this approximation is multi- 
valued, but only in its extension inside the disk, which is a non-physical part of the flow. 
At e =0 there is a transition, and the boundary of the region where the velocity potential 
becomes multivalued enters the physical space, so that the solution of Euler-Tricomi equa- 
tion cannot be considered as physically acceptable in this range of parameters. In ordinary 
viscous fluids, this multivaluedness would signal the formation of a shock wave. As said 
before, no such a thing exists in the present model of superfluid. Accordingly, the shock 
wave is replaced by a vortex, that is by adding to the flow field in the slightly supersonic 
region the flow field coming from a localized vortex close to the boundary. Taking in account 
all the nonlinear terms and the time dependent dynamics on the phase equation, we show 
that this lack of stationarity appears as the result of a saddle-node bifurcation, as suggested 



2 e is also proportional to the difference between the actual Mach number and its critical value. 
3 In this hydrodynamical limit, the microscopic length scale does not appear, and there is only 
one small parameter, e. 
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by two of us some ago |4| . We argue also that the critical velocity behaves as above the 
transonic region, being £o the healing length of NLS and R the radius of the disk; this is 
because the quantum pressure stabilizes the flow above the transonic transition for such a 
range of velocity. We claim that this saddle-node bifurcation, giving rise to time-dependent 
dynamic corresponds to vortex emission; for that, it remains to match this outer velocity 
potential, solution of the Euler-Tricomi equation with an inner solution, close to the disk 
boundary with a vortex. 

As one might expect from this rather lengthy introduction, it is a rather uneasy job to 
put together all this picture, particularly because it depends on properties of solutions the 
Euler-Tricomi equation that are not obvious. In section ||, we shall present the general 
problem in its explicit form, that is the form of the equations as well as the boundary 
conditions and the relevant facts about the hydrodynamic limit. We finally (section |HI|) 
obtain a nonlinear equation for an additional phase which can be decomposed as the usual 
Euler-Tricomi equation plus nonlinear terms that will be treated as perturbations. In section 
[TV| , we study some properties of the Euler-Tricomi equation, particularly we point out the 
multivaluedness solutions of such an equation. In section |V|, we solve the Euler-Tricomi 
equation (including the quantum pressure term) using a Fourier transform along the x— axis 
first and finding the solution in term of the Airy function. This gives a first correction to 
the critical velocity coming from the interplay between supersonic flow and regularisation 



due to the quantum pressure. Then, in section [VT|, we are able to compute solvability 
conditions around this critical value of e that lead to time dependent amplitude equation. 



These amplitude equations describe a saddle-node bifurcation. In conclusion, in section [VII] , 
we try to match this phase approach, which formally cannot describe vortices in NLS, with 
the release of vortices observed in simulation using the full nonlinear Schrodinger equation. 



II. FORMULATION OF THE PROBLEM 

The Gross-Pitaevskh equation reads, in a dimensionless form: 
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^(x, t) = -- VV + ^(x, *) |V(x, t)| 2 . (1) 

This is a partial differential equation for a complex wave function ip{x,t). This equation 
is among other properties conservative and Hamiltonian. The ground state in a periodic 
box is the homogeneous solution: ip = y/p^e~ tpot . Long wavelength and low amplitude 
perturbations propagate with the sound speed c s = y/po, on the other hand; £o ~ ^= = — 
is the only characteristic microscopic length contained in this equation. 

Writing ip = p 1 / 2 e i ^, we obtain two "real" hydrodynamical fields, p and <fi representing 
respectively the particle fluid density and the velocity potential: 

dtp=-V-(pV4>); (2) 

^ = -^W) 2 + ^vy/ 2 - P . (3) 

The first one is the density mass conservation equation, identifying by the local 
velocity v. In the second equation, the term ^T72V 2 p 1//2 , often called quantum pressure is 
negligible for large scale flows, that is for flows with a space scale much larger than the 
intrinsic microscopic length, £o- When this quantum pressure is neglected, the equation 
for cf) is the equivalent of Bernoulli's equation for a compressible fluid, with an equation of 
state for the pressure (p): p — \p 2 . For a stationary flow, d t <fi is a constant, defined by 
the conditions at infinity. Therefore, the mass density p can be computed everywhere as a 
function of v, the modulus of the velocity from this Bernoulli equation. 

^From equations @ and (|3]), we obtain the equations for the stationary flow around a 
disk (of radius R, much bigger than any intrinsic length scale £o) with a velocity at infinity 

Voo- 



V-(p(|V0|)V0) =0; (4) 

P(|V0|)=p o + i(^-|V0| 2 ) (5) 

h ■ V0 = on the disk (6) 

= v^x at infinity, (7) 
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here n is normal to the disk perimeter. 

Rewriting the equation (|4]) in the frame reference defined by the pole of the disk as the 
origin, one gets: 

d v (p(v)v)d xx (j) + p(v)d yy (f) = 0, (8) 

where x is the local coordinate tangent to the main flow, and y the orthogonal one. At 
low velocities this second order partial differential equation is elliptic. It is also possible to 
show, via an hodograph transformation, that the maximum local speed for a flow around 
an obstacle, occurs on the boundary of the obstacle. This is in some sense a nonlinear (but 
still elliptic) generalization of the min-max theorem by Riemann and Liouville. Equation 
(0) becomes hyperbolic beyond a critical velocity. This happens when d v (p(v)v) vanishes^, 
that is when the mass current takes its largest possible value for some condition at infinite. 

The above criteria (d v (p(v)v) = 0) gives, for the present model, a critical velocity v c 
such that: v 2 c = §Po + \ v 1o- When increases the property of ellipticity of equation (g) 
is broken first at the pole of the disk, leading to the nucleation of two vortices, one at each 
pole. As time goes on these vortices are convected downstream by the mean flow. These 
vortices, once released, induce as well a counterflow because of the circulation condition 
and this counterflow reduces the velocity on the surface of the disk. This brings back 
the local velocity at the pole of the disk below the critical speed, and restores there the 
ellipticity of the equation for the velocity potential. But the vortex is pulled farther and 
farther downstream, and the counter streaming effect diminishes, till the velocity at the pole 
reaches eventually the critical value, the conditions at infinity being kept constant; then new 
vortices are emitted, etc. This describes a more or less periodic release of vortices from the 
obstacle ||. We shall study in this article the process of nucleation, namely the way in 
which a vortex is emitted from the boundary when the local velocity changes slowly from 



4 This equation becomes hyperbolic also when p(v) vanishes, however, this happens for a larger 
value of the speed with the present equation of state (relation between pressure and density). 
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below critical to above critical speed. For this we developpe a phase-dynamics approach for 
the longwave asymptotics (distances larger than £ ). 



III. THE EULER TRICOMI EQUATION NEAR THE TRANSONIC REGION. 



Suppose that the local velocity at the pole of the disk, v , is close from v c \ taking this point 
the origin of the axis, the x-axis being tangential to the disk, and the j/-axis perpendicular. 
One writes the phase near the pole as = v x + yX; X being small; this gives for equation 
(|j), together with (|5|) (after an elementary redefinition of variables): 



(e + d x x)d xx x + dyyX = 



with e ~ (f o — v c )/v c . The boundary condition (|6]) becomes: 

(V V \ x^ 

vo + y<9 x x, -^dyX ) =0 at y = -— ; 

where (a, h) ■ (a', V) = aa! + bb' . Neglecting d x x in the boundary conditions gives: 



(9) 



x 



x 



d y x = —M— at y — -. 



(10) 



R * 2R 

With M = ■^ a , a constant proportional to the actual Mach number. Let us also notice that 

2 

the boundary condition arises on y = — f^, which is a parabolic approximation of the disk 
near the pole. 

Equation (^|) derives from a variational principle^ with an energy 



E = J dxdy [-~(e + ^ X ) 3 + \{d yX f 
The solution of (§) satisfying the equation and the boundary condition is 



Xo = -M-|, 



(11) 



5 This is a direct consequence from equation (Q) which is the extremum of E = ^ J dxdyp(vY 
\ I dxdy \p + i(<4-(V0) 2 ) 



S 



however, this one is not sufficient to determine the complete flow in particular the transition 
to supersonic flow, it is necessary to go up to next order. Writing \ — Xo + V 9 ; one finds: 

- (e - M-|) d xxV + d yy tp = (12) 

with 

3 

d v <p = -M^, at y = 0. (13) 



The resulting equation (|12|) for ip is Euler-TricomiQ with the boundary condition (|13|). 



The Euler-Tricomi equation may be interpreted as follows: — (e — Mj|) represents a 
generic tangential velocity profile of an ideal flow near a body, since the local main speed 
diminishes as y increases, that is as one moves far-away from the obstacle. The Mach 
number is exactly one at y = ejj|. In this equation we have neglected the nonlinear term 
dxfdxxip besides yd xx (p in (|T2"D; on the other hand the border of the obstacle has been taken 
at y — 0. (The curvature of the obstacle brings a correction to the boundary condition for 
the dominant order solution, that transforms itself to an extra term in the partial differential 
equation for the perturbation with a flat boundary <p.) This assumption is consistent to the 
following scaling in e for the coordinates 

Re'' 2 Re , Re 11 / 2 

x rsj y ~ — and in ~ - — . 14) 

M M M 3 y 1 



Once again, there exists a particular solution of equation Q12|) satisfying the boundary 
conditions, so that we will have to focus on the homogeneous solution of its equation. It 
happens that, <f> , defined as: 

,^x 3 y , r xy 3 ,^ 9 ^V 4 
(h = -M— f - eM -2- + M 



R? R? 2R A 



6 Note that here the variables are in the physical space and not the hodograph variables as in ||. 
With this procedure we have considerate directly the boundary conditions, something difficult to 
work with in the hodograph plane. 
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satisfies equation (|12|) with the flT3|) boundary conditions. 

At this point, one can iterate the linear Euler-Tricomi equation by the same procedure, 
that is by considering how <po modifies the true boundary condition and then we compute 
the next term (pi, and so on. It is then possible to find a polynomial expansion satisfying the 
full boundary condition. However, the convergence of this expansion is not obvious. This is 
not an objection in principle because we are looking only local solutions, however we have 



in mind an outer asymptotics which will mtach with a vortex kind solution (see section fVT l 



IV. SPECIAL SOLUTIONS OF THE EULER TRICOMI EQUATION 

Let us first make an "aparte" by looking at some special solutions of equation flT2|) 
that might give some idea of what arises when crossing the critical velocity, that is at the 
transonic transition. 

The roots z(x,y) of the cubic polynomial (as well as any linear combination of these 
three roots): 

, /M\ 1/3 / R\ 

z + 3 [-fri (y- e T7) z + 3x ( 15 ) 



\RJ V M) 

are exact solutions of equation (|12|). Note that for y > e||, the cubic polynomial has only 
one real root for all values of x, whereas, for y < ejj there are three real roots inside a 
semi-cubic parabola defined by \x\ < 2/3 (e — y) 3 ^ 2 . 

This multivaluation of the real roots of the cubic equation (|15D means that it is not 
possible, generally, to follow continuously a root of (|15|) along a closed path around the 
origin (see figure (|2|)). More precisely, such solution will admit a discontinuity in the region 
inside the semi-cubic. Unless one can regularize the discontinuity which arise along the 
region multivalued (in the same sense that we have to regularize the over turning of waves 
dynamics, solutions of nonlinear and non dispersive wave equations), there is no hope of 
having a stationary solution of our problem except if the discontinuity gap is 2tt. In this 
case (jump of 2ti of the phase), we will see that even if both the quantum pressure and the 
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limited transonic region have a tendency to restore the smoothness of the solution, there is 
a critical velocity above which there is no more possibility of having stationary solutions. 
The way we will treat the equation will hide the multivalued property that we pointed out 
because we will look on regular stationary solutions (which exist as well as no discontinuity 
appears), so that the discontinuities will be solved via the general time dependent nonlinear 
problem. 

A general solution arises whenever one consider a continuous superposition in the neutral 
translation mode of the above solution fll5|) : 

equivalent after a change of variable to the integral expresion in §118 of ||. This general 
expression is valid only for y > e-0 (as well as the roots of the polynomial flT5|)) where such 
a change of variable remains well defined. 

Before finish this section, let us mention, that another family of solution appears when 
one takes the Fourier transform of the x variable. Then one has that 

with $(?/) is the Airy function, the detailed analysis of this kind ok solution is elaborated 
below in a more general way. 

V. REGULARIZATION OF THE SHOCK SOLUTIONS OF THE 
EULER TRICOMI EQUATION BY THE QUANTUM PRESSURE. 

As soon as the solutions of the Euler-Tricomi equation appears to be sharper and sharper 
that the quantum pressure term in the Bernoulli equation (|3]) is no longuer negligible because 
it involves higher order derivatives. From the full Bernoulli equation one gets the value of p 
by an implicit relation (we shall consider here also the role of non stationary dependence of 
phase and density in order to capture the full short wavelength dynamics): 

P = Po + ^(^-(V0) 2 )-9 t 0+l(Ap-^), 
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reminding that = v$ ■ x + y (xo + 0o + V 9 ) an d V 9 being function of both time and position 
whereas vq, xo and 0o are independant of time. One can then estimate the value of the 
quantum pressure at the first order of perturbation, taking p = p + — (V0) 2 ) — y<9 t y? 
(it gives then the first non-zero contribution of the quantum pressure); coupling this with the 
(now) non-stationnary mass conservation equation and restoring the first nonlinear terms 
as well as the constant terms gives for the phase equation: 

- (e - M^-] d xx ip + dyyip - $d x 4(p = \d tt ip + —d tx <p + d x ipd xx (p + d x (j) d xx (p + d xx (j} d x (p + d x (j) d xx (j) . 
\ R) v l c v c 

(16) 

We have kept in this equation only the most important term of each contribution. The 
quantum pressure term, Qd x np should be multiplied by a number that we have taken to one 
by simplicity. 

The way this equation (|T6| ) is written is dictated by our method of resolution: the left 
hand side will be in fact treated as the main equation, linear and homogeneous, while the 
right hand side corresponding to perturbations which will be incorporated terms by terms. 
As said before the boundary conditions are taken homogeneous: 

dyip = at y = 0. 

Let us study the regularisation of the solution above the threshold, first by adding the 
quantum pressure then by considering cross terms involving (p and <fio; at these points the 
analysis remains linear so that we will just be able to look to homogeneous solutions of our 
problem, without solving for the amplitude; finally we will focus on the global nonlinear 
problem which allows to calculate amplitude equations and leads to the study of the time- 
dependent evolution. 

First we look at the left hand side of equation (|T6"D, taking the right hand side as zero. 
The boundary conditions are d y tp = at y = 0, so that one can expand the solution as 
a linear superposition of functions. The equation can be solved as the Euler-Tricomi one, 
using the Airy function. 
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We seek solutions of the form tp v = e x ((y), ((y) satisfying the Airy equation: 



which non-divergent solution is known to be the Airy function $(•); therefore the solution 
reads: 

A ■ e ivx § 



v 2 M 



R 



y 



R i 



2^2- 
0, 



being A a complex amplitude fixed by nonlinearities at next order. As the function Q(s) 
does not possess extrema for s > the boundary condition might be satisfied for e > 
only, otherwise A = 0. Let s n be the n th zero of $'( s ) then the only possible values for the 
wavenumber v are such that they satisfy a "quantization condition" for a given e: 



' M£ 
R 



(17) 



which have been represented on figure []. So, if e is less than a critical value, one can observe 
that, because of the quantum pressure, the homogeneous Euler-Tricomi equation has just 
the null function as solution (A = 0). This critical value (the minimum of the curve plotted 
in figure (0) is easy to evaluate: 



4x |-|Vx ^° 



4£qZ/^ with v c 



3 

Si \ 8 



M 



x 



R ^ c c V 37 \R!$J 

the critical wave number. Generically, this means that for e < e c the stationary solution 
(</>o) might describe the dynamics at the order of the Taylor expansion whereas for e > e c 
the amplitude of a solution for v — v c can expand. At this point, the amplitude cannot 
be known and has to be found as the result of the nonlinear and time dependent analysis. 
It follows now a general scheme, which consists of the evaluation of the amplitude and the 
corrections of the general solution by writing a solvability condition for their existence. 
We seek a solution of the form: 



4> = A(x) ■ e 



±iu c x 



R 



y 



3Re c 
AM 



Vi(j/)e : 



±iv c x 
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Now A(x) is a slightly varying amplitude on the horizontal variable x (the time dependance 
will be taken account in the next section) (pi being a small correction to <p depending on 
y only (the x dependance coming from the main term e tVcX ). This first correction of our 
equation ([16]) will takes in account the cross term between (p and 0o! the slow behavior of 
A(x) on x requires: 

d„A(x) 



3Re c 



4M 



A{x) 

In addition let us define the ((y), the Airy function that we use: ((y) = $((^^)^(y — 
)). The cross contribution to the nonlinear term <p x (j) xx which we have neglected until 



now, give two terms linear in ip in equation ([16|) : d x (f)od xx tp and d xx (j)od x ip. 
Therefore we obtain for the phase the following equation: 

C (fi = Wc(-( e ~ e c) - d x (f) )A(x) + iv c d xx (f) A(x) + 2zz/ c (^ - M— + d x <p )A'(x) 

+ d xx (P A'(x) - (| + M| - d x ct> )A"{x)]ay) (18) 

where Cq is the linear operator acting on the one variable function space: 

Notice that ((y) is in the kernel of Cq. The solvability condition says that the right hand side 
of (|T8|) is orthogonal to members of the kernel of Co joint operator. With the scalar product 
{fi9) — Io° f(y)9(y)dy the linear operator £ is self-adjoint, therefore ((y) belongs to its 
kernel. The solvability condition of ([R]) gives an equation for the slowly varying amplitude 
(keeping the first order in A only) A(x): 

One recognize the equation of the quantum harmonic oscillator; being the equivalent of 
the energy. It has non zero solution exists if e > e' c only. The lowest value of epsilon with 
non-zero solution satisfies: 

e c - e c = V 6e c — < e c . 

K 
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Note that this correction is small with respect to the first one. The corresponding solution 
at e = e' c reads 

(e' r -ec)x 2 2 

A(x) = Ae 8S o = 

being I the characteristic length of A(x), I = 2^f£^R/(Qe^ (note that this agrees with the 
condition of validity of this WKB approach v c l ^> 1). As for the former treatment, a non 
zero amplitude A can arise only if e > e' c ; this amplitude and its dynamics will be obtained 
by the time dependent nonlinear system for e ~ e' c . 

VI. AMPLITUDE EQUATION FOR THE SADDLE NODE BIFURCATION 

In this section we shall look at the general spatio-temporal dependence of the solution 
of ( |16D for e ~ e' c ; writing now the Ansatz for the phase ip as (A(t) = a(t) + i(3{t))\ 



(p = (a(t) sm(v c x) + (3(t) cos(z/ c x))e %*C(y) + <p\{y) cos(u c x + Q.) + ip 2 {x)((y), 

being (pi(y) and <p2{x) small corrections to (p. Imposing a solvability condition (one first 
along y, then along x), one gets the following dynamical system: 

^a"(t) - M^f3\t) = z/ c 2 (e - e' c )a(t) + <5e"^ + u 3 c {a 2 {t) - /? 2 (t)) e -^ 

1 2 2 

-(3"{t) + M^a'(t) = v 2 c {e - e' c )(3(t) + 2v c 3 a(t)/3(t) e - V (19) 

V c V c 

with 5 ~ ^/R 2 . For convenience, and because it does not change the meaning of the 
dynamics, we have taken the ratio between the different constant coming from solvability 

integration equal to one. The transcendental terms (in e~^' 2 ) are due to the interaction 

_ * 2 

between the nonlinear term and the constant term d x (pod xx (pQ with the modes sm(i/ c x)e ^ 
and cos(z/ c x)e . 

The system (|T9"|) could be written in a single complex equation for the complex amplitude 



(after an appropiate change of variable and rescaling): 
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Z tt + itoZt = fi + Z 2 . (20) 

Here \l = — (e — e' c ) 2 e^5~/4z/^ + 5e s - /^- 

The stationary solutions are Z± = ±-^/— // that is for negatives values of /z one preserves 
an odd symmetry of (p(x, y) along the x— axis, this symmetry is broken as soon as fi changes 
sign and the further evolution is more complex. 

In terms of the physical parameters the region \i < is for (e — e' c ) > — e ~ , leading 
to a new critical velocity e t , defined as: 




Notice that for small aspect ratio, the main correction for the critical velocity is in J ^ and 
is determined by e c whereas e' c and e t induce just exponentially small corrections. 

Figure shows the real roots of Z± as a function of or better as a function of 
e — et- One can identify two branches: one for low velocities (e < et); the other one for high 
velocities. As we have conducted our calculations, starting from stationary solution for low 
velocity and then adding the nonlinear dynamics step by step, this second branch for high 
velocity is irrelevant in our problem. As shown on figure (f|), for e < et there are two roots of 
the stationary equation, roots that greather for e = e 4 , whereas there is no more roots for e 
greater than e t , giving rise to a saddle-node bifurcation. The linear stability analysis of the 
stationary solutions Z± gives the following dynamic (Z(t) — Z± + z(t) , and \z(t)\ \Z±\): 

z tt + iuzt = ±2^— \iz. 

Therefore, one notes that for: 

one branch is stable (Z_, representing by A_ on figure (^)) and the other is always unstable 
as in ordinary saddle-node bifurcations. However, because of to the oscillatory term iuzt, 
when approaching the bifurcation, both solutions are stable. This happens in the range of 
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values of e such a that 




< e < e t 



The resulting phase diagram of the bifurcation has been plotted on figure @. We argue 
that this saddle-node bifurcation gives a consistant scenario of the vortex nucleation seen 
in figure (p. They appear as a consequence of the disappearence of stationnary solutions 
on the phase dynamic approch. The vortex diminishes the local velocity so that the flow 
comes back to a description valid for e < et and when the vortex is far enough, because of 
the advection due to the mean velocity, e pass again through the transition and we got a 
periodic vortex-nucleation process. The saddle-node bifurcation appears in our case to be 
richer because close to the transition both stationnary solutions are stable. This result is in 
a good agreement with numerical solution by Huepe and Brachet [|J . 

VII. BREAKDOWN OF THE PHASE DESCRIPTION, THE APPEARANCE OF 
VORTEX MOTION, CONCLUSIONS AND FURTHER MISCELLANEOUS 

We have in the former section reviewed the phase description, that explains the dis- 
appearence of any stationary solution for the flow problem. This is the first step toward 
the nucleation of a vortex. Briefly, one can retain that this periodic behavior appears as 
a saddle-node bifurcation where two branches of stationnary solutions collapse giving rise 
for larger velocities to non stationnary solution at all. The numerous critical parameters 
we have mentionned might be simplified if one see that as R 3> £o> we have et ~ e' c ~ e c . 
This allows to claim that the critical velocity for vortex nucleation v v reads at first order 
correction of the critical speed: 



where v c is the Landau critical velocity. 

Finally, it remains to match the vortex nucleation as a process being part of the same 
evolution of the dynamical saddle-node bifurcation. For simplicity we will consider the 
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amplitude equation (|20|) where we have dropped the first order time derivative, a term much 
smaller than the other ones as one approaches the time when nucleation occurs (see the 
scaling below): 

Z{t) = Z 2 (t) + ^i 

where /i is related to e — et. When \x is negative, there is two real stationnary solutions (Z = 
the stable, z = \J—[L the unstable) corresponding to a±. The saddle-node bifurcation 
is crossed when /i becomes positive. This can be studied by taking \i = t (by rescaling, no 
multiplying factor is needed) and then we obtain the first Painleve transcendent^] II: 

Z = Z 2 + t (21) 

A convenient change of variable, for t < is: 

z = ^W(T); T=~(-t) 5 / 4 

5 

and gives the following equation: 

W 4 W o 

w + = W 2 - 1. 

T 25 T 2 



7 One may note that taking \x = t means to take e = et in equation (|l6|). This gives 



et - d xx (p + dyyip - \d u <p - —d tx f = 0, 



which, after a change of variables of the form: r\ = et — My/R and ( = Mv^t/R + ey/vQ, and 
neglecting the linear derivative in time as in (H), leads to: 



-V9xx<p + (M 2 /R 2 - e 2 /vl) (d m <p - d a <p) = 0. 



A kind of Euler-Tricomi equation, interesting by itself. One note that, if the rate of acceleration 
at infinity e is larger than Mvq/R, the nucleation process is caused by a dynamical instability not 
by a sonic transition as is described in this article. 
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It typically gives the two stationnary solutions W = ±1 for large T and it can be solved 
in terms of Weierstrass function for large T ||. Equation ([H]) is known for giving finite- 
time singularities which shows that the slowly varying approximation of the phase dynamics 
breaks at some point. In addition, one can notice that Z is related to the velocity of 
the fluid in the moving frame (of velocity Vq). Then, one expect to relate the amplitude 
Z to a order parameter which parametrizes a continuous family of solutions of the full 
nonlinear Schodinger equation (fj). Jones and Roberts [JTo|] found the kind of solution that 



we are interested, in two and three spatial dimensions, consisting in axisymmetric solitary 
structures. The solution proposed is 

<P = -Ux + m ) J 22 
x 2 + (1 — U z )y z 

m being a constant and U the free parameter characterizing the whole family of solutions. 

Whenever the parameter (a velocity) U of these solutions is close to the unity (the sound 
speed) one identifies the relative speed U — 1 (via a change of variable) with the true order 
parameter of the saddle-node bifurcation. 

On the other hand as U goes to zero one matches the solution ( |2"2"D with a pair of moving 
vortices, one with negative topological charge located in (x = 0, y = a) while the other is an 
image one inside the disk (x = 0, y ~ — a)f| 



Ux , ((y + a)VT=TP\ f(y- a )VT^TP^ 
arctan — arctan 



Vl — U 2 \ x J \ x 



y/l^TP x 2 + (l-U 2 )y 2 ' 

One relates very easily the distance a in ( p3|) with the speed U in ( p2]) by U = j-' ^ ne 
Hemholtz law of motion for point-like vortices. As a diminishes, a vortex pair appears. 
So in the dynamical context a satisfy formally the same equation than the quantity z of 
(PT|) or a of flTP|), when the time dependant solution evolves (after crossing the saddle-node 



^Note that at first order the dependence on the radii R of the disk is not relevant to this solution. 
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bifurcation), a vortex appears as a increases, the second vortex of the pair being formally 
inside the disk in other to preserve the boundary condition. 

Finally, it is of a general interest to note that the shallow water equations have exactly 
the same shape that the set (|2|,|3|) in the case of potential fluid motion (see 0). In fact the 
shallow water equations are 

d t h = -V- (for); 
d t w + v • Vv = —gSJh. 

After imposing a vortexless flow: rot v = 0, i.e. v = V0, and neglecting the quantum 
pressure term in (||), one identifies the height of fluid h with the superfluid density p, and 
the potential fluid velocity with the phase of the condensate wave function. It is possible to 
get a short scale term in the shallow water equations by adding a capillary term, therefore 
it exist a complete analogy with the analysis developed in extenso in this article. Perhaps a 
transonic transition with the predicted scenario is observable in mercury where the kinematic 
viscosity is very low. However, the final state will be different because the deep significance 
of the phase of the wave function ip does not extend the existence of quantized circulation 
vortices for the case of surface waves, where as it is well known the circulation takes any 
value. 
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FIGURES 




FIG. 1. Numerical simulation of the nonlinear Schrodinger equation for a bidimensional flow 
around half a disk; the velocity at infinity is Voo = 0.442 and we have taken dx = 0.125 the mesh 
grid and the radius of the disk is R = 7.5. a) & b) respectively the modulus and the phase of the 
wave function at t = 20 time unit of NLS. The densite and the phase go up from bright to dark 
color. One can see the low density around the top of the obstacle, due to a Bernoulli effect, c) <fe 
d) same functions at t = 50.6. A low density structure is advected by the flow (at right of the top 
of the disk). One can see a phase discontinuity and the tip of it where the phase is not defined this 
is the signature of a topological defect, that is a quantized vortex. 
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FIG. 2. Shape of z(x,y), for e = around the origin. 
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FIG. 4. Stationary amplitude of the phase equation for ^ = 0.1. The part for e > e' c is not 
valid in our approximations whereas we have dashed the unstable solutions. 
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FIG. 5. Stability diagramm of the stationary solutions a± in the lu complex plane (the stability 
is studied via the expression w(t) = a-e tujt , a being an amplitude. The small dashed line corresponds 
to the real and imaginary axes whereas the large dashed line indicates the stability evolution of 
a + as e increases untill e< where the stationnarity disappears as the collapse of a + and a_ . The 
thick line represents the same evolution for a_ . 
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picture vortex 

FIG. 6. Position of the vortices for the unstable stationnary solution. 
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